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It is well known that loss of information about a system, for some observer, leads to an increase 
in entropy as perceived by this observer. We use this to propose an alternative approach to de- 
coherence in quantum field theory in which the machinery of renormalisation can systematically 
■ be implemented: neglecting observationally inaccessible correlators will give rise to an increase in 

D ' entropy of the system. As an example we calculate the entropy of a general Gaussian state and, 

assuming the observer's ability to probe this information experimentally, we also calculate the cor- 
£f~^ ■ rection to the Gaussian entropy for two specific non-Gaussian states. 

PACS numbers: 03.65.Yz, 03.70.+k, 03.67.-a, 98.80.-k 

or 

^ , I. INTRODUCTION 

The most natural way of defining the entropy of a quantum system is to use the von Neumann entropy: 

> 

Os\ S vN = -Tr[pln{p)}, (1) 

where p denotes the density operator which in the Schrodinger picture satisfies the von Neumann equation: 
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g: i4 t p = [H ,p], (2) 



where H is the Hamiltonian. Since quantum mechanics and quantum field theory are unitary theories, the von 
Neumann entropy is conserved .albeit in general not zero. 

The decoherence program [1(3 usually assumes the existence of some environment that is barely observable to some 
^ ■ observer. This allows us to construct a reduced density operator p roc i characterising the system S, which is obtained 
by tracing over the unobservable environmental degrees of freedom E: /5 rc d = Tr^;[/5]. This is a non- unitary process 
which consequently generates entropy. This operation is justified by the observer's inability to see the environmental 
degrees of freedom (or better: access the information stored in the ES-correlators) . Yet, tracing out some degrees of 
freedom results in complications. The simple looking von Neumann equation ^ is replaced by an equation for the 
reduced density operator. This "master equation" 0, E| can be solved for only in very simple situations, which has 
hampered progress in decoherence studies in the context of interacting quantum field theories. In particular, we are 
not aware of any known solution of the master equation that would include perturbative corrections and implement 
the program of renormalisation. Also, one has no control of the error in a calculation in the perturbative sense induced 
by neglecting non-Gaussian corrections. 

Here we propose a decoherence program that can be implemented in a field theoretical setting and that also allows 
us to incorporate renormalisation procedures. The idea is very simple, and we present it for a real scalar field (j)( x )- 
A generalisation to other types of fields, e.g. gauge and fermionic fields, should be quite straightforward. 

The density matrix p(t) contains all information about the (possibly mixed) state a quantum system is in. From 
the density matrix one can, in principle, calculate various correlators. For example, the three Gaussian correlators 
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are given by: 

Wx)4>(y)) = Tr[p(t)4>(x)$(ff)}=F(x,t-,y,t')\ t=tl (3a) 
= Tr[mHx)Hy)} = d t d t 'F(x,t-,y,t%= t , (3b) 

i({0(f),7r(y)}) = ^'IV[p(t){^) 1 7r(^}]=^F(x ) i;f,0|t=t' J (3c) 

where n denotes the momentum field conjugate to <p, and where the curly brackets denote the anti-commutator as 
usual. Note that these three Gaussian correlators can all be determined from the statistical propagator F(x, t; y, t') 
in the Heisenberg picture as usual in quantum field theory. The fourth Gaussian correlator is trivial as it is restricted 
by the canonical commutation relations. There is an infinite number of higher order, non-Gaussian correlators that 
characterise a system, where one can think of e.g.: 

■ ■ ■ $(x n )) = Tr\p(t)$(Sx) ■ ■ ■ fan)] (4a) 
(n(x 1 )---Tr(x n )) = Tr[p(t)7r(afi) • • • n(x n )] , (4b) 

where n > 3. Moreover, one can also imagine a higher order correlation function consisting of a mixture of both </>'s 
and 7r's. Finally, one could think of other non-Gaussian correlators obtained by further differentiating the statistical 
propagator: 

d?d?F(x, t; y, t% =t , = d?dFTr[p(t )4>(x, t)4>(y, 011*=*' , (4c) 

where n + m > 3. For free theories, all the non-Gaussian correlators either vanish or can be expressed in terms of the 
Gaussian correlators. 

Let us turn our attention to interacting field theories and present two simple arguments why, generally speaking, the 
Gaussian correlators dominate over the non-Gaussian correlators. In interacting quantum field theories, the Gaussian 
correlators all stem from tree-level physics, whereas non-Gaussian correlators arc generated by the interactions. Rely- 
ing on a pcrturbativc treatment of these interactions, we can safely expect that all of these higher order, non-Gaussian 
correlators above are suppressed as they generically are proportional to a non-zero power of the (perturbatively small) 
interaction coefficient. A second, more heuristic, argument invokes the (classical) central limit theorem, according to 
which systems with many stochastic mutually independent degrees of freedom tend to become normally distributed, 
and hence nearly Gaussian. While this theorem is derived for classical systems and for independent stochastic vari- 
ables, we believe that it will in a certain quantum disguise also apply to weakly interacting quantum systems with 
many degrees of freedom. 

Hence, we expect that to a good approximation, many of the relevant properties of quantum systems are encoded 
in the Gaussian part p g of the density operator p and the non-Gaussian correlators are thus typically much more 
difficult to access 1 . The crucial point is that all information about a general Gaussian density matrix and thus about 
the Gaussian part of a state, is stored only in the three equal time correlators ([3]). All non-Gaussian correlators 
like ((4]) contain information about the correlations between the system and environment as well as the environment's 
correlations induced by higher order interactions. When projected on the field amplitude basis in the Schrodinger 
picture, the Gaussian part of a density matrix reads: 

p s [ct>,(t>';t} = (<t>\p s (t)\(f>')=Me X p -| dlxd^ (^(a: )i(£,y;t)^(y ) + f (x)B(f, - 2#)C(!:, y;t)^'(y)) 

(5) 

where A,B,C and N are functions of the equal time field correlators (for details see sections [Vj. Of course, the 
von Neumann equation ^ does not hold any more for p g , and thus the entropy conservation law is violated by p g . 
Indeed, neglecting the information stored in observationally inaccessible correlators will give rise to an increase of the 
entropy 2 . 



1 As a simple example, consider the temperature correlations induced by scalar field perturbations from inflation: while the amplitude 
of scalar field fluctuations is given by Tr[p(t)ij>ff (x,t) 2 ] ~ GffH 2 ~ 10~ 12 , where </>jv denotes the Newtonian potential, where Gjv is 
Newton's constant and H is the Hubble rate during inflation, its quantum corrections are expected to be of the order ~ (G^H 2 ) 2 ~ 10 — 24 
or smaller. 

2 In passing we note that we shall not address entropy generated due to the loss of unitarity caused by some observer's inability to access 
the information stored in all of the space-time. Prominent examples of this research field include black hole and de Sitter space entropy 
and there is an extensive literature on the subject. In [f| truncation of the infinite hierarchy of Green's functions was proposed as an 
operational definition of the entropy of cosmological perturbations. Calzetta and Hu [f| coined the expression "correlation entropy" and 
proved an H-theorem for the correlation entropy in a quantum mechanical model. 
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The purpose of the present work is to calculate the entropy of a system taking these considerations into account. 
We distinguish two cases. In the simplest case we assume that our observer can only measure Gaussian correlators 
and his apparatus is insensitive to all non-Gaussian correlators. We will consider this situation in sections ITTl and ITTll 
In quantum field theories, the problem seems to be that the evolution equations of the three Gaussian correlators ([3]) 
do not close as soon as perturbative corrections are included. The reason is that off-shell physics becomes important. 
A way out of this impasse, and in our opinion the way out, is to solve for the dynamics of the unequal time statistical 
correlator, i.e. the statistical propagator, from which one then extracts the three relevant correlators @. 

In the second case, we assume that our observer is sensitive to specific types of non-Gaussianities present in 
our theory, apart from the leading order Gaussian correlators of course, either because the observer has developed a 
sensitive measurement device or the non-Gaussianities happen to be large enough. We will consider this case in section 
IIVI The sections above contain, for pedagogical reasons, quantum mechanical derivations rather than quantum field 
theoretical calculations. We will generalise our results to quantum field theory in section fVl 

This paper aims only at developing the necessary machinery to discuss entropy generation in an interacting quantum 
field theory in an out-of-equilibrium setting, see e.g. 0-dj]. We discuss the application of these ideas elsewhere [13, 13 
in great detail, where we consider a scalar field <f> coupled to a second scalar field via the interaction Lagrangian density 
£i nt = — ^h^ 2 {x)(f){x) . The results we develop here are applicable to a much more wider class of field theoretical 
models, including gauge and fcrmionic fields. 



II. GAUSSIAN ENTROPY FROM THE WIGNER FUNCTION 



Entropy has of course been widely studied over the years [19h24| . In this section we will derive an entropy formula 
based on a naive probabilistic interpretation of the Wigncr function. This characterises much of the 20th cen tury 
efforts to connect the Wigner function, Boltzmann's equations, Boltzmann's H-theorem and macroscopic entropy [251 ] . 
Let us begin by considering the Hamiltonian of a time dependent harmonic oscillator, which is the one-particle "free" 
Hamiltonian: 

H(p, x) = £- + -muj 2 x 2 + H S , H s = xj . (6) 
2m 2 

In general, the oscillator's mass and frequency and the source current can all depend on time: m = m(t), uj = io(t) 
and j = j(t), respectively. H s is the source Hamiltonian and j = j(t) the corresponding current. This Hamiltonian 
is relevant both for many condensed matter systems that can be realised in laboratories and in an early Universe 
setting. For the latter, the time dependence of the parameters is introduced for example by the Universe's expansion 
or by phase transitions. Alternatively, this Hamiltonian can represent for example the simplest model of a laser, where 
x = x(t) is the photon amplitude, where, if m — > 1, oj(t) — > E(t)/H is the photon's frequency and where j = j(t) is 
a c harg e current that coherently pumps energy into the system, transforming the vacuum state into a coherent state 

IE if. 

The time dependence of the parameters changes the system's energy, but simultaneously does not change the 
Gaussian nature of a state if such a state is imposed initially, i.e.: the evolution implied by the Hamiltonian 
transforms an initial Gaussian density matrix into another Gaussian density matrix. Moreover, as we will come to 
appreciate, one has to add a coupling to the environment or add a self-interaction to generate entropy. The Lagrangian 
L = px — H follows as usual from the Hamiltonian. Note that the source term can be removed from this Lagrangian 
by a simple coordinate shift: 

x(t) ->z(t) = x(t) -x (t), (7) 
upon which the Lagrangian reduces to a quadratic form: 

L = 2 ?n ^ 2 — —muj 2 z 2 + Lo(t) + — [mzxo] ■ (8) 

Here, Lo(t) is a z-independent function of time and the last term is a boundary term which does not contribute to 
the equation of motion if xo(t) is a function that obeys: 

^- [mx ] + mu} 2 x + j = . (9) 
at 



A general Gaussian density matrix centered at the new origin remains a Gaussian density matrix centered at the origin 
under the evolution of a quadratic Hamiltonian and moreover, which is what the analysis above shows, linear source 
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terms in the Hamiltonian will not alter the Gaussian nature of the state. We can thus consider Gaussian density 
matrices centered at the origin, whose time evolution is governed by an Hamiltonian of the form ^ with j(t) — > 0. 
When written in position space representation, the single particle density operator of a quantum mechanical Gaussian 
state centered at the origin is of the form: 



Ps(t)= dx d y\x}p g (x,y,t)(y\, (10a) 

J — oo J — oo 

where: 

p s (x,y;t) = jV(i)exp [-a(t)x 2 - b(t)y 2 + 2c(t)xy] , (10b) 

and where a = a(t), b — b(t) and c = c(t) are determined from the von Neumann equation ([2]). The subscript g 
denotes "Gaussian" . Moreover, from p' g = p g it follows that b* = a and c* = c. When c = one recovers a pure state 
with vanishing entropy. However when c ^ the density matrix is mixed and entangled and it cannot be written 
in the simple diagonal form, p g (x,y;t) = ^*(y,t)^(x,t), where ^(x,t) — VJ7 cxp(—ax 2 ). When c > the density 
matrix tends to get diagonal in the x — y direction, however when c < the diagonalisation occurs in the x + y 
direction in which case the entropy is not defined. The Heisenberg uncertainty relation is in trouble in this case as 
we will sec in equation (|30[) . Hence it is natural to assume that c > 0. We keep the discussion phcnomcnological and 
do not discuss the physical origin of c ^ 0. The normalisation Af is obtained from requiring Tr[/3 g ] = 1: 

/°° f°° I 7T 
dx(x\p g \x) = dxp £ (x,x;t)=J\f x — r = l, (11) 
-oo J-oo V A a R _ c ) 

from which we conclude: 



Af=f-^^ t (12) 

provided that c < or where clr = 3? [a]. The equations that the functions a(t), b(t) and c(t) of the density matrix (I10b|) 
obey can easily be obtained from the von Neumann equation (J5|), which in this Gaussian case in the amplitude basis 
reads: 

h 2 1 
iHd t p g (x,y;t) = -— (d 2 x ~ d 2 ) p g (x, y; t) + -mu 2 {x 2 - y 2 )p e {x,y-t) . (13) 

If we insert equation (|10b|) in the equation above we find: 

— - = — aia R (14a) 
dt m 

dC AH /HUN 

— = — oic 14b 

dt m 

In (TV) - — ai (14c) 

dt m 

- + (14d) 

where we defined a\ — 3 [a]. Note that Af cx y/oR — c is consistent with equation (|14c|) as it should 3 . 
The Wigncr function is defined as a Wigncr transform of the density matrix p g (x,y;t): 

/oo 
dre-^/ H p g (q + r/2,q-r/2;t), (15) 
-oo 

where we defined the average and relative coordinates, q — (x + y)/2 and r — x — y, respectively. Hence, a Wigncr 
transform can be thought of as an ordinary Fourier transform with respect to the relative coordinate r of the density 
matrix. The inverse Wigner transform determines p g (x,y;t) in terms of W(g,p;t): 

p g (x,y;t) = J"" ^Lert'-vVtyvfant), (16) 



3 Indeed, by making use of equations l|27|l one can show that equations l|14j l are consistent with (and equivalent to) the Hamilton equations 
for the correlators 11471 . 
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where 2ttH is the standard phase space measure. We obtain the Gaussian Wigner function by performing the inte- 
gral ([15]): 

W(q,p;t) =M{t)exp[-a w (t)q 2 - (3 w (t){p + q Pc (t)) 2 ] , (17) 

where: 



a w = 2(aR, — c) (18a) 

= ^R-M (18b) 
p c = 2ha\ (18c) 



and p c denotes a "classical" momentum, which makes the Wigner function non-diagonal. Note that the Wigner 
function (|17|) is normalised to unity, in the sense that: 

!>(,,p; t ) = l. ^ 
, znn 

Recalling definition (JTSJ), equation (|13l) transforms to: 

f ft + -ft - ™ 2 ?3 P ) V%, p; t) = , (20) 
where we made use of the following relations: 

(^-^)/>g(z,2/;t) = 2d q a rPg {q + r/2,q-r/2;t) (21a) 
re Wi = _ lfidpe tpr/h ^lb) 

Upon noting that p/m = x = d p H and —miu 2 x — p = —d x H, we see that equation (|20p is nothing but the Liouville 
equation for the Boltzmann's distribution function f(x,p;t), 

d t f(x,p;t)+xd x f(x,p;t)+pd p f(x, P ;t) = 0. (22) 

A probabilistic interpretation of the Wigner function thus seems plausible by the following identification: 

W(q,p;t)^f(x,p;t). (23) 

Of course, the equations for yV(q,p;t) and f(x,p;t) are identical only in a free harmonic oscillator theory and, 
when interactions are included, differences arise between the von Neumann equation for the density operator p g (t) 
(or the corresponding equation for W(q,p; t)) and the Boltzmann equation for f(x,p;t). This indeed makes the 
identification p3]) less rigorous but it remains a useful picture. The correspondence (|23|) has been widely used in the 
literature to develop a formalism relevant for example for the physics of heavy ion collisions [TH, [28| or for baryogenesis 
HQ. 

In the distant 1969 Cahill and Glauber [27j have proposed to approximate the density operator in the coherent 
state basis \a) by its diagonal form: p g w p cs (a) = (a\p g \a), which they interpret as a quasi-probability distribution. 
Since coherent states harbour many properties of classical particles, one can argue that this is reasonable. Recall 
that coherent states are defined by a\a) = a\a), where a = \J mw/2{x + ip/{mw)] denotes the annihilation operator 
for the oscillator ([6]). The density matrix p C s(u) represents an early example of a quasi-probability distribution for a 
quantum system. Needless to say, by this line of reasoning one can get only approximate answers for the expectation 
value (A) = Tv[p g A] of an operator A and likewise for the entropy (which is just the expectation value of — ln(p g )). 
Another common way of approximating expectation values of operators is to give a probabilistic interpretation to 
the Wigner function, along the same line of reasoning as argued above in equation (|23[) for a free theory. In this 
approach [13, HH the entropy is approximated by: 

S w S w = -Tr[W(q,p; t) ln(W(q,p; t))] , (24) 
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where the trace Tr — > J dpdq/[2irh] should be interpreted as an integration over the phase space volume as in 
equation fTS]). We now insert {T7| into ([H} to find: 

S W = ilnf^)+l-ln(2) 
2 \a R -cJ 

= In f ^ + 1 , (25) 



(26) 



where we defined: 

A 2 = «R + c 
br-c 



We evaluated the phase space Gaussian integrals by shifting the momentum integration to p' = p — p c (the Jacobian 
of this shift equals unity). At this point, it is useful to evaluate a few quantum mechanical expectation values. We 
can extract the following correlators from our Gaussian state: 



i 2 ) = Tr[p g x 2 ] = / dx(x\ Ps x 2 \x) = (27a) 

J-oo 4 (aR - c ) 

\{*,P}) = - K, ^ , (27b) 
2 / 2(a R - c) 

(f) = H 2 Z. c2 . (27c) 

or - c 

Here, we made use of (x\p\ip) = — ihd x (x\ip) . Moreover, one can easily verify that ([x,p]) = (xp — px) = ih as it should. 
This enables us to find the equivalent inverse relations: 

ai -^rn^r (28a) 

A 2 + l 

flR = w ( 28b ) 

A 2 - 1 



and to express A(t) in equation in terms of the correlators (f2"T|) : 



(i 2 )(p 2 >-(i{i,p})^ 



(29) 



For the moment this is just a definition, but as we will come to discuss, the physical meaning of hA(t) is the phase 
space area in units of h occupied by a (Gaussian) state centered at the origin. For a pure state A(i) = 1, while for a 
mixed state A(i) > 1. Hence, we can define the uncertainty relation for a general Gaussian state: 

(hA(t.)) 2 ,„ 9w , s / 1 rn ,\2 h 2 

{ I" =(£ 2 )(f)-(^{i,p}) > T - (30) 

If c is negative and |c| < or, we can in principle have A < 1, which can be seen from equation (|26[) . It would 
be interesting to investigate whether any physical experiment could be proposed where a violation of Hcisenberg's 
uncertainty relation occurs, even if it were for only a very brief moment in time (such that it would still hold on 
average) . 

Furthermore, it is natural to define the statistical particle number density on our phase space, in the quantum field 
theoretical sense, as: 

, , Am - i 

n{t) = -ii , (31) 

in terms of which the Wigner entropy (|25j) can be rewritten as: 

S w (t)= In (n(t) + ^\ +1. (32) 
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These results in essence agree with the entropy per mode of scalar cosmological perturbations obtained in [3l|. 
Shortly we shall see that the Wigner entropy (|32[) is a good approximation for the von Neumann entropy if n 3> 1, or 
cquivalently when h/S. 3> 1. 

Let us try to put these results briefly in historical context. The fundamental question that is left unanswered 
in [3~lj is what the dynamical mechanism for the growth of A is. Instead, the authors argue that an observer that 
measures cosmological perturbations can be approximated by a coherent state (which is a good classical basis) , thus 
explaining the growth in A. Two mechanisms were considered in [32| : a projection on the particle number basis and 
on the coherent state basis. In both cases, information of the phase of the squeezed Gaussian state is lost, which thus 
generates entropy. In the latter approach formula (|29|) is reproduced, while in the former, the entropy acquires the 
form of the statistical entropy of Bose particles S = (n + 1) ln(n + 1) — n ln(n) with n — > sinh 2 (r) the particle number 
associated with the squeezed state, and r the squeezing parameter of the state. While these early articles, and many 
subsequent ones [HI, [34j] , attempted to explain the origin of entropy of cosmological perturbations by a (non- unitary) 
process of coarse graining (associated with particular properties of the observer's measuring apparatus), a genuine 
dynamical mechanism of entropy generation was lacking. 

This has changed recently [EL ItR I35l - l40j . now dynamical mechanisms for entropy generation have been proposed. 
Rather than discussing various mechanisms in detail, let us just point out that the notion of how entropy should 
be defined and generated has evolved: recently, it has been advocated that the relevant entropy to consider is 
the Gaussian entropy, although the approaches differ in the actual origin of the entropy increase. For example, 
Prokopec and Rigopoulos suggested in [37J that tracing over the unobserved isocurvature perturbations leads to 
entropy generation in the adiabatic mode. Similarly, tracing over the unobserved gravitational waves should lead to 
entropy generation in the adiabatic mode (4l| . Kiefer, Lohmar, Polarski and Starobinsky [H, Hsj as well as Campo 
and Parentani @, among others, suggested that, if super- Hubble modes take the role of the system, the origin of 
the entropy of cosmological perturbations should be associated with the stochastic noise generated by the sub- Hubble 
ultraviolet modes representing the environment. Both Giraud and Serrau [42j and we [13, 5H H] advocate that 
neglecting obscrvationally inaccessible non-Gaussianities is a prominent mechanism to generate entropy. 

So far, we have discussed how the Gaussian entropy in Wigner space is defined and qualitatively discussed various 
mechanisms to generate entropy dynamically. Let us now, to finish this section, discuss how one can sensibly think 
of a decohered state in Wigner space and discuss the physical meaning of A(f). Despite the fact that the Wigner 
entropy (|25|) only approximates the von Neumann entropy, it offers a useful, intuitive way of depicting decohered 
states irrespective of the precise underlying mechanism by which such a state has decohered in the first place. Let 
us, for this reason, discuss this phenomenologically. The Wigner function (|T7|) is a 2-dimcnsional Gaussian function 
whose width in the p-direction need a priori not necessarily be the same as in the g-direction. Hence, in general it 
can be squeezed in some direction. This is illustrated in figure [T] The la or 2<r cross-section of the phase space 
distribution has the shape of an ellipse, which we also visualised in figure [T] If the cross-section is a circle centered at 
the origin, the state is either pure or mixed and if, moreover, its area is unity A(t) = 1 it is a pure vacuum state. If 
the cross-section is again a circle but displaced from the origin, we are dealing with a generalised coherent-squeezed 
state (whose shape is not described by equation (f2"5)) for obvious reasons). If the cross-section is an ellipse, it is a 
squeezed state as mentioned before. The ellipse is parametrised by: 

a w (t)q 2 + p w {t){p + qp c (t)) 2 =0, (33) 

where we have made use of equation (|I7[) and where $ is some constant that determines the height at which we slice 
the Wigner function. The area of this ellipse is now given by 4 : 

A(t) = d Z = torhA® , (34) 
^a w (t)/3 w (t) 

where we have used equations l|18p and (|26p . Clearly, the area of the ellipse in Wigner space is determined straight- 
forwardly from the phase space area A(f). Note that if = 1/tt we have A = A. 

Figures [2] and [3] visualise the decoherence process of a squeezed state. Initially, in figure[2j the state is squeezed and 
has unity phase space area A (to) = 1, such that its entropy vanishes. Now we switch on our favourite decoherence 
mechanism: the state interacts with some environment and environmental degrees of freedom are traced over or 
non-Gaussian correlators generated by the interaction with the environment are neglected. Consequently, the area 
in phase space increases A(t) > 1, such that S(t) > 0. Hence, the ellipse in Wigner space grows which we depict 
in figure [3J An important feature of any decoherence process reveals itself: for a highly squeezed state, given the 



4 Note that the area of ellipse defined by the equation Ax 2 + Bxy + Cy 2 = 1 is A- = 2n/\/ iAC — B 2 . 
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Figure 1: General shape of the Wigner function, rep- 
resenting the 2- dimensional Gaussian phase space dis- 
tribution of a certain state. The la or 2a cross-section 
is an ellipse. For the current quantum mechanical dis- 
cussion, we of course let 7r — > p and (f> q. 




n 




Figure 2: Squeezed state before decoherence. For ex- 
ample, the area of the ellipse of a vacuum state would 
be unity A = 1. For the current quantum mechanical 
discussion, we of course let -k — > p and <f> — s> q, also in 
figure 



Figure 3: Squeezed state after decoherence. Clearly, 
the accessible phase space has increased and one can 
say that knowledge about the system has been lost. The 
entropy has increased compared to figure [J| 



knowledge of q after some measurement, the value p can take in a subsequent measurement is constraint in a very 
narrow interval. Interestingly, after the state has decohered, and given the same measurement of q, the range of values 
p can take has increased. Indeed one can say that knowledge about p has been lost compared to the squeezed state 
before decoherence, hence entropy has been generated. Alternatively, we can say that given some measurement of 
q, the value p can take is, after decoherence, drawn from a classical stochastic distribution, which indeed coincides 
with the familiar idea that decohered quantum systems should behave as uncorrelated stochastic systems. Indeed, n 
roughly counts the number of patches in phase space that behave independently and arc uncorrelated. 

With this intuitive notion of decohered state in mind, one can understand that, for highly squeezed states, the 
position space basis (which generalises to the field amplitude basis in the quantum field theoretical case) is the 
pointer basis. Let us start out with a highly squeezed state, i.e.: take for example a pure state during inflation that 
rapidly squeezes during the Universe's expansion. The Hamiltonian of the state is then dominated by the potential 
term, and the kinetic term only contributes little. In thermal equilibrium, the system will minimise its free energy 
F = H — TS, where T is the temperature of the heat bath. Now, if we switch on some decoherence mechanism due to 
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interaction with the environment, the entropy will increase mainly due to momentum increase, i.e.: (p 2 ) will increase. 
Indeed, increasing (p 2 ) will hardly affect the Hamiltonian, whereas it will significantly affect the entropy. Hence, the 
variable x is robust during the process of decoherence, in the sense that (x 2 ) will hardly change such that it qualifies 
as a proper pointer basis. Note however, that x is only a pointer basis in the statistical sense, such that there is a 
well defined probability distribution function from which a measurement can be drawn. 



III. GAUSSIAN ENTROPY FROM THE REPLICA TRICK 

The entropy of a general Gaussian state has been derived by Sohma, Holevo and Hirota 0, |4|| by making use of 
Glauber's P representation for the density operator. Here we present an alternative derivation for the von Neumann 
entropy (Q} of a general Gaussian state by making use of the replica trick (see e.g.: [47]-[4!|). Of course, there are 
other analogous methods one can use 5 . Notice first that: 

ln(p g ) = lhn^-i. (35) 

Here, e is a positive integer which is analytically extended to zero. Hence, in order to calculate the entropy ([T]) we 
need to evaluate Tr [/S g +C ] ■ Using equation (jTUJ) , we thus have: 

/OO /"OO 
dx • • • / dx e exp [-2a R (x 2 H h x\) + 2c {xqx-i + x x x 2 H V x e x )] 
-OO J — OO 

/OO pOO 
dx --- / dx e exp [—2/3 {(x Q — axi) 2 + (x\ — ax 2 ) 2 H h (x £ - ax Q ) 2 }] 
-OO J —OO 

d2/exp[-2/3y 2 ] 



1+e 

' —oo 

1+e 



= |l-a 1+ r 1 (^J ■ (36) 

Here, J = |1 — a 1+e | _1 is the Jacobian of the transformation from x t to = Xi — ax i+ i, where i = {0, 1, • • • , e — 1}, 
and y e = x e — axo. A comparison of the first and the second line in (|36p tells us that: 



c ca^ 



In the limit c — > one must recover a pure diagonal density operator and, consequently, the shift a should vanish, 
singling out: 



8 = 13- . (38b) 

in equation (|37p as the physical choice. Based on (a R — c)/(3 = (1 — a) 2 , we can rewrite equation (|3!)|) in a simpler 
form: 

Recall that normalisability of the density operator requires or > c, implying that < a < 1, such that the absolute 
value in (|39[) can be dropped. This means that Tr [/5g +e ] is analytic in a (for positive integers e): 

T r [Pl +e ] = > (40) 



5 We thank Thco Ruijgrok for his useful comments. 
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such that one can analytically extend to complex e. According to the replica method, the entropy is then obtained 
by taking the limit e — > 0. To perform this procedure we expand (|40[) to linear order in e and we find the following 
expression for the von Neumann entropy (JlJ: 

5 vN (t) = -ln(l-a)--^ln(a). (41) 
1 — a 

Using the expression for a in equation (|38a[) and the definition of the phase space area in equations (j2"6")l and (j2T))) , we 
find: 

such that we can express the entropy (|4"Tj) solely in terms of A: 



Since < a(t) < 1, we have 1 < A(t) < oo. The von Neumann entropy (|4"3"]l vanishes when A = 1 (or c = 0), defining a 
pure state, while a strictly positive entropy 5 v n > implies a mixed state A > 1. Clearly, the Gaussian von Neumann 
entropy is completely determined by the three Gaussian correlators characterising the state (|29|) . Relations (|4"3]) 
and (|41|) suggest the following definition of the particle number n(t): 

„ (()s ^) = ^|Z1, (44, 
1 — a(i) z 

where < n < oo, in terms of which the entropy (|4"Tj). (|4"3")l becomes: 

SVnW = (1 + ln(l + n) - nln(n) , (45) 

which is the well known result from statistical physics for the entropy of n Bose particles per (quantum) state j50| . 
Note that (|45[) is a convex function of n, such that S^rt!) + S(ri2) > S(ni + n 2 ), which is another desirable property for 
the entropy. The particle number n defined in (|44[) should be interpreted as the number of independent (uncorrelatcd) 
regions in the phase space of a single (one particle) quantum state. 

We will show next that A is a special function, as it is conserved by the evolution equations resulting from the 
Hamiltonian ©: 

'x = d p H = — (46a) 

TO 

p = -d x H = -mw 2 x , (46b) 
where we set j — > 0. These equations imply for the correlators: 

= (47a) 

J t {l&P}) = -^ 2 (x 2 ) + l(p 2 ) (47b) 

±([x,fi) = (47c) 

±<f) = -2TO W 2 (i{i,p}). (47d) 

Recall that these equations are, in the light of equation (|28p. equivalent to (fT4")) . The third equation is automatically 
satisfied by the commutation relations, [x,p] = ih. One can combine the other three equations in (|47p to show that: 



d h 2 A 2 d 



dt 4 dt 



{x 2 ){p 2 )-(\{x,p})' 



= (48) 



This implies that A (and in fact any function of A) is conserved by the Hamiltonian evolution 

Of course, when interactions arc included, the conservation law (|48|) becomes approximate. In quantum field theory 
interactions are typically cubic or quartic in the fields and will thus induce non-Gaussianities in the density matrix. 
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Figure 4'- Left: The von Neumann entropy (solid line) and its Wigner entropy approximation (dashed line) 
for a Gaussian state as a function of the phase space area A. Right: The difference between the Wigner and 
the von Neumann entropy. The Wigner entropy approaches the exact von Neumann entropy rapidly as A 
grows. 



In particular, for a quantum mechanical A</> 4 model, this has been investigated by Calzetta and Hu [61 l5lj. They 
derive an H-thcorcm for A in the case of a quantum mechanical 0(N) model. 

Finally, it is interesting to compare the von Neumann entropy (|43|) and the Wigner entropy (|25p . as shown in figured 
This reveals that the Wigner entropy represents a large phase space (semiclassical) limit of the von Neumann entropy, 
which should not come as a surprise. Indeed, the Wigner entropy can be used to accurately represent the entropy in 
systems which develop large correlators (occupation numbers), such as cosmological perturbations [3l|. The Wigner 
entropy fails however to take account of quantum correlations present in the density matrix and Wigner function, 
which are properly taken into account in the exact expression for the entropy of a general Gaussian state (1431) . 

IV. NON-GAUSSIAN ENTROPY: TWO EXAMPLES 

An important question is how to generalise the von Neumann entropy of a Gaussian state to include non-Gaussian 
corrections. Non-Gaussianities can be important either because they are made large on purpose by specially preparing 
the state, or simply because they are accessible in measurements due to the fact that the observer's measurement 
device is very accurate. 

While we are far from building a general theory of how non-Gaussian correlations affect the notion of phase space 
volume, statistical particle number and thus entropy, we shall present two examples in this section: we consider a 
system in a state which can be represented by an admixture of a Gaussian density matrix and some small non-Gaussian 
contributions. Secondly, we consider small quartic corrections to the density matrix, contributing to the kurtosis of 
the ground state. 

Qualitatively, we expect that the correction to the Gaussian entropy due to non-Gaussianities present in a theory is 
small, except for systems whose entropy is (almost) zero. This can most easily be appreciated in the Wigner approach 
to entropy. Consider for example a system with a large entropy, for example a squeezed state whose phase space area 
in Wigner space has increased significantly due to developing a A » 1. Non-Gaussianities can deform the area in 
phase space of the state (by squeezing, pushing or stretching the Wigner function) roughly by a factor of order unity. 
Consequently, the effect on the entropy will be relatively small. This can be understood by realising that a lot of 
information is contained in unequal time correlators. Of course, such an argument ceases to be true for a pure state, 
whose quantum properties are very pronounced. 

A. Admixture of the Ground State and First Excited State 

Let us consider a density matrix of the form: 

p 2 (x, y;t)=M 2 (l + Cox + Co V + Ci^ 2 + CiV + ^xy) exp (-ax 2 - a*y 2 + 2cxy) . (49) 
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Such a state typically appears when one uses a laser to excite an harmonic oscillator in its ground state. The system 
can then be described by an admixture of the ground state and a small contribution to the first exited state 6 . Note 
that the process of pumping energy in a simple harmonic oscillator in its ground state could lead to creating a state 
that is not pure (c ^ 0). The measuring apparatus is assumed to be sensitive to the admixture of the two states. We 
shall assume that non-Gaussianity is small, in the sense that all Q parameters (i = 0,1,2) are small, such that we 
content ourselves with performing the analysis up to linear order in Q. 

Notice first that to this order the terms in parentheses in equation (|49p can be approximated by an exponential: 

1 + Coa: + CoV + (ix 2 + CiV + %2xy ~ exp (( x + Qv + (ix 2 + Qy 2 + 2C, 2 xy) , (50) 

implying that the first two terms (,qx + Qy induce an entropy conserving shift in x and y (as in coherent states). 
Of course, this ceases to be true at quadratic and higher orders in On the other hand, the £i and £2 terms do 
change the entropy even at linear order. At higher order, £1 and £2 induce kurtosis whose effect on the entropy we 
will discuss shortly. 

Before we turn our attention to calculating the entropy, let us first show that a density matrix (|49[) with £ 
indeed generates non-trivial higher order correlators. Normalising the trace to unity yields: 

M = V^K + C+i(ClH + C 2 )- 

An interesting higher order correlator to consider is for example the connected part of the four-point correlator: 

(i 4 )con = (i 4 )-3((i 2 )) 2 . (52) 

If the connected part of the four-point correlator is non-zero, the state is said to have kurtosis. Using the Gaussian 
density matrix in equation (|10p one can easily show that the connected part of the four-point correlator vanishes as 
it should: for free theories higher order correlators either vanish or can be expressed in terms of Gaussian correlators. 
Using however the non-Gaussian density matrix in equation (|49p with Co ~^ and the normalisation constant (|51[) wc 
find: 

/a4\ _ 3(Clfl +C2) 2 3(Cli? + C2) 2 , n , (/ . /roN, 

{X con - -—, ^7 ^ --TT-, rj" + U {{Qir) , {(,2) ) ■ (53) 

16(a R - c) 2 (a R - c + ±(& R + C2)) 2 16(a R . - c) 4 

Clearly, genuine non-Gaussian correlators are generated and this state has kurtosis at quadratic order. 

Let us now calculate the entropy by making use of equation ([5H)l . As said, it is precisely d and (2 that affect the 
entropy. To see that notice further that their effect can be captured by a shift in a and c as follows: 

a — > a = a — £1 (54a) 
c -> c = c + C2 , (54b) 

such that the corresponding von Neumann entropy simply reads (cf. equation (05])): 

S 2 = (1 + n)ln(l + n) - nln(n) , (55) 

where fi = n + Sn and where as usual n = (A — l)/2. Now Sn cx , £2 is the correction to the Gaussian entropy we 
are about to calculate. Comparing with equations (|57|) and (|4"4")l we have: 



Recall that the wave function of the pure one particle state of the simple harmonic oscillator J6J (with j = 0) is: 



i[>l(x) = (x\ipi) 



4/mo;\ 3 l 1 / 4 



x exp I 



211 



with energy Ei = frw (l + |). The problem at hand reduces to this state upon identifying art mu/(2h), Oi — > 0, c — » 0, ^0 - ► and 
Ci -> 0. 
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from which we immediately find, to linear order in £m and (2: 



6a = ?(±^l(Jl^ + l + * 2 . (5?a) 



1 + a \ or — c 2a aR — c 
6a n(l + n) f ( 1R 1 + 2n + 2n 2 C2 \ 

(1-a) 2 l + 2n Ur-c 2n(l + n) a R - c/ ' 1 j 

where we used osr/c = (1 + a 2 )/(2a). Finally, <552 = { m [(l + implies that: 

S 2 (t) = S g (t) + 6S 2 (t) , (58a) 

where 6S2 is an entropy shift given by: 

= ^ f-^- + 1 9 + ;; + 2 f C2 ) m (I±^ , ( 58b) 

w l + 2n V a R- c 2n(l + n) a R -cJ \ n J ' 

where S s = (1 + n) ln(l + n) — n ln(n) is the von Neumann entropy of a Gaussian state and n is the statistical particle 
number associated with the Gaussian part of the state as before. 

Several comments are in order. Firstly, the result (|58|) implies that to linear order 3[Ci] docs not change the 
entropy. In fact, 5[£i] induces (dc)squcczing of the state, which can be appreciated from equation (|18cj) . Secondly, 
6S2 is positive (negative) whenever and £2 are positive (negative), irrespective of n > 0. Thirdly, we have rescaled 
Cir and £2 by (or — c) to get a dimensionless quantity. This rescaling is natural, since (Zr — c measures the width of 
the state. Finally, for large n the formula (|58[) gives a meaningful answer, since: 

Um ^4^^, (59) 
2 or-c 



is finite. On the other hand, when n — >• 0, we encounter a logarithmic divergence: 



lim <5S 2 = ( + l-^—] In (- ) , 



n-s-0 \ OR — C 2 Or 

indicating a mild (logarithmic) breakdown of the linear expansion. Notice however that also in that limit the for- 
mula (|55j) remains applicable. 



B. Kurtosis 



In the former example non-Gaussianities are generated by adding a quadratic polynomial to the prefactor. One 
could also introduce non-Gaussianities by adding higher order, i.e.: cubic or quartic, powers in the exponential. Cubic 
corrections generate skewness and since they contribute to the entropy at quadratic and higher orders in the skewness 
parameter, we shall now focus on the quartic corrections to the density matrix. Quartic corrections to a density 
matrix affect the kurtosis of a state and we parametrise this as follows: 

P4,(x,y;t) =W 4 cxp (-ax 2 - a*y 2 + 2cxy) exp (770a; 4 + ?/q?/ 4 + r^x 3 ?/ + r^xy 3 + f] 2 x 2 y 2 ) = —^p g (x,y;t)p ng (x,y;t), 

(61) 

where p g is the Gaussian density matrix (|10p as before, with the normalisation A/" g given in equation (|12p . and where 
finally p ng = exp(?7oiE 4 + VoU^ + Vi x3 y + Vi x y 3 + V2X 2 y 2 )- As in the previous example, we shall consider only linear 
corrections in r)i (i = 0, 1, 2) to the entropy. To this order the non-Gaussian part of (|61[) can also be written as: 

6p m = png-i = m* 4 + Vay A + mx 3 y + vlxy 3 + m^ 2 y 2 + £>{vf) ■ (62) 

We can find the normalisation constant after some simple algebra by making use of ([62]) : 

JV 4 =JV f 

To gain some intuitive understanding of what our density matrix looks like, let us consider figure [5J Here, we show 
the Wigner transform of equation (|61[) using (1621) , with ryi = = 972. Clearly, the quartic corrections change the peak 



4 [2(o R - c)] 2 



(63) 



14 





Figure 5: Wigner transform of the density matrix described by equations \61\) and $62\) , with i)o / only. 
We used 770 > (left) and 770 < (right). The peak structure of our state in phase space, that was so simple 
for Gaussian density matrices, is much more complicated. Some regions in phase space have a negative value 
of the Wigner function, indicating a breakdown of the use of a Wigner function as a convenient measure of 
probability. 



structure of our state in phase space. Moreover, some regions in phase space now have a negative Wigner function. 
This nicely illustrates the limitations of using the Wigner function as a probability density on our phase space as soon 
as non-Gaussianitics (due to interactions) arc included, i.e.: equation (|23[) holds only approximately in this case. 
In order to calculate the entropy, we need Tr[/54 +e ] as before (cf. equation (|36p ). To linear order in Sp ng we have: 



Tr[p^ 



(64) 



The first term in the parentheses is just the Gaussian result (|36p while the latter term can be evaluated by making 
use of the formulae in Appendix [A] Expanding to linear order in e, the equation above yields: 



= 1 



ln(l - 



1 



■ ln(a) 



J Vo + Vo 
4 [2(or - c)] 2 



(1 



4a 



ln(a) 



1 



3 m +r,l /l 



-In(a)-1 )+--, 



172 



4l+a + a 2 



ln(a) - l) 



4[2(a R -c)] 2 Vl-a v ' ) 4 [2(a R - c)] 2 V3 1 
As expected, we see that the entropy naturally splits into a Gaussian and a non-Gaussian contribution: 

Si[t) = S s (t) + SSi(t) , 



(65) 



(66a) 



where S g is the Gaussian entropy as before, and where: 



Vo + Vo ( 3n(l + n) 



[2(a R -c)] 2 V 



In 



1 



[2(o R - 



1 + 2n \ 
1 + 3n + 3n 2 



3 m+jrt 

4 [2(a R -c)] 5 



(1 + 2n) In 



l + 2n 



hi 



(66b) 



This is the main result of this section and intuitive in the following sense: just as in the first non-Gaussian example 
above, we see that positive kurtosis parameters {f/oi?.7 fjijj, ^2} > 0, tend to make the effective state's width (o R — c) -1 / 2 
larger, which increases the area in phase space the state occupies, which in turn increases the entropy. If however 
{VoR-, Villi V2} < 0, the phase space area shrinks, which in turn decreases the entropy. The statements above hold for 
any statistical particle number 0<n< 00, with the exception of n — > 0, where, just as in the case studied above, 
a weak logarithmic divergence occurs when r/m ^ or 772 7^ 0. Notice finally that 3[t7o] and Ss[t)i] again do not 
participate in entropy generation, but rather contribute to the squeezing of the state. 

Kurtosis and skewness in quantum mechanics, studied in this section, occur also in interacting quantum field theories 
which we discuss next. 
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V. ENTROPY IN SCALAR FIELD THEORY 

The quantum mechanical expressions for the Gaussian and non-Gaussian entropies that we have developed for 
pedagogical reasons in sections [Til Hill and IIVI for a single particle density matrix can be generalised to field theory. 
We firstly need to consider correlators in quantum field theory however. 



A. Equal Time Correlators in Scalar Field Theory 

Let us now proceed analogous to equation (|10[) and write the density matrix operator for our system in the field 
amplitude basis in Schrodinger's picture (see e.g. HH): 



where: 



p s = Jv<pJ 2?0'|0>Pg[0,'/'V](0'h 
p g [<j),^;t] = N exp(-0 T '• A- <p - <j)' 7 ' ■ B-<j>' + 2</> T '• C*-</) , 



(67) 



(68) 



where T>4> = Yisev d4>(x) and \<p) = Ylsev \ < t > {^))- A few words on the notation first. We shall consider <fi = <j>(x) 
as a vector whose components are labelled by x. Moreover, A(x, y, t) can be viewed as a matrix, such that A ■ <p is 
a vector again, where a ■ denotes matrix multiplication which, for the case at hand, is nothing but an integral over 
D — 1 dimensional space. Hence, quantities like T • A ■ </> are scalars and involve two integrals over space. Note that 
at this point we do not assume that A is homogeneous, i.e.: A(x,y,t) ^ A(x — y,t), but rather keep any possible 
off-diagonal terms for generality. 

Our density matrix (|68| is hermitian, such that we have A*(x,y,t) = B(x,y,t) and C*(x,y,t) = C(x,y;t), where 
we also used A(x,y,t) = A(y,x,t) and C(x,y,i) = C(y,x;i). The normalisation Af of p[<f>,<p';t] can be determined 
from the standard requirement Tr[/5] = 1: 



dct 



2(A R - C) 



1/2 



(69) 



where we note that the hermitian part A^ of A is real and symmetric such that A^ = Ar. 

Just as in the quantum mechanical case, we need to calculate the three non-trivial Gaussian correlators which 
completely characterise the properties of our Gaussian state. In order to calculate these correlators, we need to have 
an expression for the statistical propagator as in equation ([3]) which, given some initial density matrix p(to), is in the 
Heisenberg picture defined by: 



Let us begin by calculating: 



F(x, t; y, t') = Tr p(t )j>{x, t)j>(y, t') 



0{$)4>{y)) - F{x,y;t) = Tr p(t)$(x)$(y) 



(70) 



(71) 



where we have made use of the Heisenberg evolution equation for operators. It is convenient to add a source current 
j{x,t) to the density matrix (|67[) . such that p becomes: 



4^, <f/;t] = Wexp (-<P T - A ■ </> - <P' T ■ tf-tf + 2<j) T ■ C-t/f + f 

in terms of which equation (JTTJ can be rewritten as: 
S 6 



>+r-4t 



(72) 



F(x, y- 1) 



Sj{x) Sj{y) 

_6 S_ 

Sj{x) 5j(y) 
1 



j=j'=o 

T>4>Afexp (-24> T -(A K -C) 



exp i(j+/) T .(4,-C)- 1 .( j +/) 



(A R -C)- 1 (x,y;t) 



3=3' =0 

(73a) 
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We also need the other correlators: 



\{{m,Ky)}) = d v F{x,t;y,t')\ t=v = ±Tr{p g (t){j>(x),n(y)}} = ~{A K - C)" 1 ■ A^x^t) 
i({7r(f) j7 r(y)}) = d t d t ,F(x,t;y,t% =t , = ^[pS){HZ),nm 



(73b) 
(73c) 



A^-iAn-Cy^A+A-iAn- Cy 1 ■ A^ - C '-(A K - Cy 1 ■ C 



-ih-Ar5[(j)' — (/)]. As a check one can verify that: 



Ti(p g (t)[<l>(S)^(y)}) = ihS^ix - y ) . 
Combining the equations above we find: 

,. , 1 



(x , y; t) . 



(74) 



A'(x,y;t) 



4 

(Ax-Cy 1 ■(A R + C)(x,y;t) 



(75) 



This is the desired field theoretic generalisation of the Gaussian invariant A 2 in equation P^]) . Notice that the result 
above applies for general non-diagonal Gaussian density matrices. This is an important quantity because, just as in 
the quantum mechanical case, A will be the conserved quantity under any quadratic Hamiltonian evolution. Since 
the von Neumann entropy is also conserved in this case, it is natural to expect that SVn = SVnIA]. 

Of course, if we are interested in problems in which the hamiltonian density is only time dependent, one can make 
use of spatial translation invariance of the correlators, such that the equal time statistical correlator is homogeneous: 
F(x, y; t) — > F(x — y, t). In this case it is beneficial to Fourier transform according to: 



F(k,t) 



l {x - y)F{x ~ y,t)e- li - {S -y ] 



(76) 



such that equation (|75l) becomes local in momentum space and reduces to the result known in the literature: 



A 2 (M) = « 



F(k,t,t)d t d t 'F{k,t,t') - (d v F{k,t,t') 



r, 2 

(A R + C)(k\t) 



(A K -C)(k,t) 

This representation is particularly useful in problems with spatial translational symmetry, such as cosmology |3lL |4 



(77) 



B. Entropy of a Gaussian state in Scalar Field Theory 

Let us now discuss the von Neumann entropy ([1} of the Gaussian density matrix (|67|) by using the replica 
method (f3"S")) . One can proceed analogous to section UTT1 Some subtleties arise however as we deal with a system 
with infinite degrees of freedom. For this reason, we nevertheless include an outline of the proof in appendix [B] The 
entropy for a quantum system that can be described by a Gaussian density matrix is given by: 



S, 



vN 



Tr 



A + I 



In 



A 



In 



A -I 



(78) 



We denote the identity matrix by I = S D 1 (x — y). As in the quantum mechanical case, we can define the generalised 
statistical particle number density correlator n = n{x, y, t) as: 



I 



(79) 



in terms of which the entropy (|78p reads: 



S vN = Tr [(n + I) ■ In (n + 1) - n • In (n)] 



(80) 
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In the limit when ||n|j 3> 1, in the sense that for the diagonaliscd matrix its diagonal entries of interest are large, 
we can expand the logarithms in to get for the entropy, SVn ~ Tr[ln(?i) + I + ©(n^ 1 )], which nearly coincides 
with the Wigner entropy Sw m held theory, cf. equation ([32|) . Equations ([80jl and (|78|) represent the von Neumann 
entropy of a general Gaussian state in scalar field theory and are the main result of this section. 
In the homogeneous limit, we can again Fourier transform and equation (|78[) reduces to: 



S y 



V 



d 3 k 

(2tt)3 



A(k,i) + 1 /A(k,i) + 1\ A(k,i)-1 /A(k,t)-1 



In 



In 



(81) 



where the volume factor arises because of the trace. The limit A(k, t) ^> 1 basically agrees with |3l| . 



C. Entropy of a Non-Gaussian state in Scalar Field Theory 



Analogous to the one particle non-Gaussian entropy discussed in section ITVl let us generalise that result to the field 
theoretical case. The non-Gaussian density matrix p 2 in equation (|67|) generalises to: 



p 2 = Jv^J 2W)p 2 [0,0';i](0'|, 



where: 



P2 [<P,<P';t] = N 2 [l + Co • </> + C f • </>' + 4> T <i-<P + cjy' T + 4> T ' • Ca 



(82a) 



6' + • C| • 0l po[0,0';i] (82b) 



p [(/),(j)';t] = exp 



■C-c 



C T ><f> 



(82c) 



We could of course perform a similar shift in the exponent as we have done before in equation (|54[) in which case 
the resulting equation for the entropy would formally be exact, i.e: we do not assume yet that the non-Gaussian 
contributions are small. If we want to expand around the Gaussian result however, we can only perform the integral 
if we assume that all correlators, including the non-Gaussian ones, are homogeneous, i.e.: they are only a function of 
the difference of their coordinates A(x, y, t) = A(x — y,t). The Fourier transform of equation (|52")l is: 

(83) 



P2^A'-A = n^kcxp [-A(k,i)|0(k)| 2 - A*(k,t)|0'(k)| 2 + 2C(M)0*(k)0'(k)] +0((f) , 



where we have set Co = as before, as it will not induce an entropy shift. Also, the product over the momenta is only 
over half of the Fourier space. Moreover, we have absorbed the small non-Gaussian contributions in the functions A 
and C as in the quantum mechanical case: 



A(k,t) -+ A(k,i)=A(k,t)-Ci(M) 
C(k,t) C , (k,t) = C(k,i) + C 2 (k,t), 

Now we can read off the result for the entropy in Fourier space in equation J5TJ) as: 



S y 



V 



d 3 k 

(2tt)3 



A(k,i) 



l ln ^A(k,i) 



1 



A(k,t) 



i ln /A(k,t) 



2 V 2 j 2 V 2 

Assuming that Q(k,t) <^ 1 we can derive the change in entropy to linear order in Q. The result is 

S 2 (t)=S g (t) + SS 2 (t), 
where is the Gaussian contribution and 6S2 is an entropy shift given by: 



SS 2 



V 



d 3 k 

(2^)3 



Cifl(k,t) 



l + 2n(k,t) + 2n 2 (k,t) 



C 2 (k,f) 



A R (k,t) - C(k,t) 2n(k,t)(l + n(k,t)) A R (k,t) - C(k, t) 
y n(k,t)(l + n(k,t)) fl + n(k,t) 



1 + 2ra(k,i) 



n(k,t) 



(84a) 
(84b) 



(85) 



(86a) 



(86b) 



where, of course, n(k, t) = (A(k, t) — l)/2. This result is the field theoretic generalisation of the quantum mechanical 
entropy (|58|) . It can be applied to mildly non-Gaussian states in field theory, which, apart from being mixed, also 
contain small one particle contributions. 

The field theoretical generalisation of the second example presented in section HVBI is hard to solve for, even in the 
homogeneous case, as it is non-local in Fourier space. 
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D. Example: A Scalar Field with a Changing Mass 

As a simple illustration of the ideas presented above, let us investigate the effect of a changing mass on the Gaussian 
entropy for a scalar field. This is maybe not a very exciting example, as no entropy is generated of course in free 
theories, it nevertheless provides an intuitive way of how the Wigner function can be used. Let us consider the action 
of a free scalar field: 

S[</>] = J d 4 x l-ld^Wd^Wrr ~ \rnl{t)4> 2 {x)} , (87) 

where as usual rj^ = diag(— 1, 1, 1, 1) is the Minkowski metric, and where we consider the following behaviour of the 
mass m,p(t) of the scalar field, mediated for example by some other Higgs-likc scalar field: 

ml(t) = (A + Stanh(pt)) . (88) 
The equation of motion following from (|87[) reads: 

(d 2 -d 2 + m 2 (t))cb(x)=0 1 (89) 
Let us quantise our fields in D-dimcnsions by making use of creation and annihilation operators: 

^ = / (kMt>** + • (so) 

The annihilation operator acts on the vacuum as usual a%\0) = 0. We impose the following commutation relations: 
[dg, Og] = (2ir) D ~ 1 6 D ~ 1 (k — k'). The mode functions 4>k(t) of 4>(x), defined by relation thus obey: 

(a* + w 2 (t))fo(f) = o, (91) 

where k = \\k\\ and ui 2 (t) = k 2 + m^(t). Using equation (|7Tj) with p(to) = |0)(0|, we see that the mode functions 
determine the statistical propagator completely: 

F^(k,t,t') = ^{0fc(f )#(*) + &(*)#(*')} • (92) 

The solution of (|91[) which behaves as a positive frequency mode in the asymptotic past, i.e.: limt_>_oo = 
exp [— iuJi n t] /^/2w; n , can be expressed in terms of Gauss' hypergeometric function 2-F1 (see [53l |54|): 



0*(*) = -/==«P 



—iuj + t — i — log{2 cosh(pi)} 
P 



2 F 1 ( l + t^ ) ,^;l-»^;i{l+tailh(/>i)}) , 

P P P 2 



where we defined w- m = (k 2 + A — B)z , uj out = (k 2 + A + B)? and w± = (w ou t ± Wj n )/2. Having the mode functions 
at our disposal, we can find the rather cumbersome expressions for the exact statistical propagator. The statistical 
propagator in turn fixes the phase space area through equation (|77|) . yielding A&(t) = 1 such that: 

S k (t)=0. (94) 

We thus conclude that a changing mass does not change the entropy for a free scalar field. Let us now examine the 
same process in Wigner space. Of course, we will reach a similar conclusion (|94j) but Wigner space is much more 
suited to visualise the process neatly. The statistical propagator is also the essential building block for the Wigner 
function which can be appreciated from generalising equations p8|) and (fT8|) to: 



aw(M) = 2F„(k,t,t) (95a) 

MM) = 2F tf£ t] (95b) 

M k (t) = (95d) 
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Figure 6: Squeezing of the Gaussian phase space in the 
Wigner representation due to a changing mass. We 
used k/p — 0.2 and m/p changes from to 4 where the 
mass changes most rapidly at tp = 5. For early times 
tp — 0, the Gaussian phase space is a perfect circle 
(blue). Already for tp — 2 a little squeezing is visible 
(red) and at late times tp = 9 this is manifest (black). 
The coordinates p r and q r have been rescaled for di- 
mensional reasons: p r = p\Ju(t) and q r — q/ ^Juj(t). 



Figure 7: Rotation of the Gaussian phase space in 
the Wigner representation due to a mass that has 
changed. We used the same parameters as in figure fj>| 
The squeezed phase space ellipses are plotted at times 
tp = 9.5 (green), tp = 9.6 (blue), tp = 9.7 (red) and 
tp — 9.8 (black). The rotation is clearly visible and its 
direction is indicated by the arrow. This can be under- 
stood from realising that an arbitrary squeezed state can 
be interpreted as a superposition of ordinary coherent 
states and the knowledge that the latter show the same 
rotating behaviour. 



We can thus plot cross-sections of the 2-dimcnsional Gaussian phase space distribution in the Wigner representation. 
In figure [5] we depict the squeezing of the Gaussian phase space due to the changing mass. At early times, the state 
is still in an unsqueezed vacuum. At tp — 2 some squeezing is visible, whereas at late times tp = 9 this is manifest. 
Despite the effect of the changing mass on the accessible phase space in Wigner space, the area of the ellipse remains 
constant throughout the whole process as expected from equation (|M| . At late times, when the mass has settled to 
its final constant value m/p — 4, the squeezed ellipse rotates in Wigner space, which we depict in figure [71 but it is 
not squeezed any further. The rotation can be anticipated by the intuitive notion that an arbitrary squeezed state 
can be thought of as a superposition of coherent states (just imagine that the ellipse is replaced by a number of circles 
displaced from the origin). Coherent states that are displaced from the origin rotate in time. 



VI. CONCLUSION 



We have formally developed a novel approach to decoherence in quantum field theory: neglecting observationally 
inaccessible correlators will give rise to an increase in the entropy of the system. This is inspired by realising that 
correlators are measured in quantum field theories and that higher order n-point functions are usually pcrturbativcly 
suppressed. An important advantage of this approach is that the procedure of rcnormalisation can systematically be 
implemented in this framework. 

We have shown how knowledge about the correlators of a system affects the notion of the entropy associated to that 
state in two cases: we firstly calculated the entropy of a general Gaussian state. This entropy can be expressed purely 
in terms of the three equal time correlators characterising the Gaussian state. Moreover, all three correlators can 
be obtained from the statistical propagator which opens the possibility to study quantum corrections on the entropy 
in an interacting quantum field theory in a systematic manner. Secondly, we calculated the entropy for two specific 
types of non-Gaussian states. Firstly, we assumed that the state of the system can be described by an admixture of 
the ground state and a small contribution of the first excited state. This yields a small correction to the Gaussian 
entropy. Secondly, in the quantum mechanical case, we calculated an expression for the entropy in case our observer 



20 



could probe a specific type of kurtosis of the ground state. This also changes the entropy that the observes associates 
to the state. 

We also outlined the use and limitations of the Wigner function, the Wigner transform of the density matrix of 
a system. In particular, we have shown that A, the fundamental quantity constructed from various correlators that 
fixes the entropy, indeed coincides with the phase space area in Wigner space. 

This is a rather phcnomcnological discussion, connecting the notion of entropy, correlators and phase space area 
in quantum field theory. Although we have outlined various mechanisms how entropy could be gen erated, we have 
in the present paper not applied our ideas to concrete systems. We refer the reader to e.g. fl7l. Tl8l. |42| for specific 
examples of how entropy is generated in interacting quantum field theories in an out-of-cquilibrium setting. 
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Appendix A: Useful Integrals 

Here we quote some integrals that are used in section [IVJ The three integrals are: 



V(0 



Me) 



dx ■ 



dx e [r) x A t + ^Xq) exp [-2a R (xo H h^) + 2c(x xi H h x e x )] 

(l-a 2+2£ ) 2 



4(2/3)^ (l-a 2 ) 2 (l- a i+ e )5 

/GO 
dx e (r/iXgXo + r]lx e xl) exp [—2a-&(xl H h x 2 e ) + 2c(x xi H + # e #o)] 
-co 

3 ^ a(l + a- 1+g )(l-a 2+26 ) 
- W 1+r 'i) 4(2/3) ^ (l- a 2 ) 2 (1 _ a i +e) 4 

/CO 
dx e {rf 2 x 2 e xl) exp [-2a R (xQ H h x 2 ) + 2c(x xi H h x e x )] 
-oo 



dx 



^27 



1 7T^" (l+a 2 )2(l+ a -2+2 £ )_ Q ,4( 1+Q ,-6+2, :)+6a 



(2/3Y 



(l-a 2 )2(l- a i+^3 



(Al) 



(A2) 



(A3) 



where e > is an integer. The first non-trivial check of these integrals is to set e = 0, in which case their sum gives 
the (inverse) normalisation constant (|63l) . 



/,o(0) + /,i(0) + /„2(0) 



~3r)o+Vo +Vi +V* +V2 



2(a R - c) 4 [2(o R - c)] 2 



= N7 1 - AC 1 



Moreover, we encourage the reader to check the integrals for some other values of e. By performing direct integrations 
we have checked the integrals (|A1HA3|) for e = 1,2,3,4,5, and they all give the correct answer. Since the expres- 
sions (| All I A3[) are analytic in e when < a < 1, they can be uniquely analytically extended to complex e's in that 
range of a, based on which we calculate the von Neumann entropy in section HVl For the entropy calculation we need 
to expand the integrals (|A1HA3|) around e = up to linear order in e. From equation (|64p we see that it is convenient 
to multiply the integrals with Af} +e = [2(a R - c)/7r] (1+£)/2 , 



AC+^oW = 7 
A/- g 1+ %r(e) = - A 



3 ?7o + Vo 



4 [2(a R - c)] 2 

3 Vi+Vl 

4 [2(a R - c)] 2 

3 rj 2 

4 [2(a R - c)] 2 



1 + e 
1 + e 
1 + e 



, /-, x a (5 + at) , , , 

ln(l - a) + — ^ ln(a 

1 — a z 

1 + 3a + 2a 2 



ln(l - a) + 



1-a 2 



ln(a) 



i / 1 \ 14 + 7a + 7a 2 , , ' 
ln(l - a) + — — 2 In(a) 



(A4) 
(A5) 
(A6) 



where we used a R /c = (l + a 2 )/(2a), 2/3 = c/a = 2(a R — c)/(l — a) 2 , and we rescaled r}i by the Gaussian width of 
the state squared, (a R — c) 2 to write it in natural dimensionlcss units. 
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Appendix B: Von Neumann Entropy for a Gaussian Field Theory 



Here we outline how to calculate the von Neumann entropy in field theory for a quantum state that can be described 
by a general Gaussian density matrix. Keeping equation (|35[) in mind, we realise that, as in the quantum mechanical 
case, we have to evaluate: 



Tr [p 1+e ] = N 1+t / V(j) ■ ■ • / V<t> t exp 



T>4> e exp 



T>4> e exp 



= N 1+t J Vfr 



2(4% -A R - 0o + ■■■ + £■ A n ■ 4> e ) 
2($'-i R I-& + ■■■ + $■ A K I-4> e ) 



(Bl) 



where we note that Ar and C are real and symmetric matrices and we assume that they can both be diagonalised 
by the same orthogonal matrix O. We transformed to the diagonal field coordinates by setting e.g. </>o = O(j>o, such 
that all matrices diagonalise, e.g.: O -1 • An ■ O = Ar(x, t)I. Moreover, we assumed that C T = C. Secondly, we 
transformed variables to <pi = 4>i — al^i+i for i = {0, 1, • • • , e — 1} and <p t = (f> e — al</>0j analogous to the quantum 
mechanical case. The Jacobian of this change of variables is now given by: 



where we write down the matrix R for clarity as: 



Idetf.R- 1 ] I =det (I-al)" 



(B2) 



R = 







al ■ • • \ 
I -al ■■■ 
I ••• 



V -al 

The diagonal matrices a(x,t)I and /3(x,t)I generalise to: 

/ 

a(x, t)l = 



(B3) 



I / 



Ar(x, t) 



P{x,t)I 



( A R (x,t) S 

C(x,t) \\C(x,t) / 
\ 

C(x,t) i 
2a(x,t) ' 



1 I 



where we note that Ar — C — (3(1 — a) 2 . Equation (|B1[) can now easily be evaluated: 



Tr [p 1+t ] = dot 



(!-") 
1 - a e 



e+l 



Expanding around e = 0we get: 



-Tr 



ln(l - a) + 



1 - a 



ln(a) I 



(B4) 
(B5) 

(B6) 
(B7) 



where we made use of, det[I + eA] = 1 + eTr[A] + C(e 2 ). The phase space invariant A 2 (x,y,t) in equation (|75|) can 
also be cast in diagonal form: 



A»(j;t)i=^Mi 

A R (cc,t) - C(x,t) 



(B8) 
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We can thus derive a(x,t) = (A(x,t) — l)/(A(x, t) + 1). The resulting equation for the entropy can be rotated back 
to the original non-diagonal field coordinates by making use of the unitary matrix O: 



S v n = Tr 



A + I / A 
— r — • In 



In 



A -I 



(B9) 



which is nothing but equation ((751) . 
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